General features of the energy landscape in Lennard-Jones like model liquids 
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Features of the energy landscape sampled by supercooled liquids are numerically analyzed for 
several Lennard-Jones like model systems. The properties of quasisaddles (minima of the square 
gradient of potential energy W = jVFj 2 ), are shown to have a direct relationship with the dy- 
namical behavior, confirming that the quasisaddle order extrapolates to zero at the mode-coupling 
temperature Tmct- The same result is obtained either analyzing all the minima of W or the saddles 
(absolute minima of W), supporting the conjectured similarity between quasisaddles and saddles, as 
far as the temperature dependence of the properties influencing the slow dynamics is concerned. We 
find evidence of universality in the shape of the landscape: plots for different systems superimpose 
into master curves, once energies and temperatures are scaled by Tmct- This allows to establish 
a quantitative relationship between Tmct and potential energy barriers for LJ-like systems, and 
suggests a possible generalization to different model liquids. 



I. INTRODUCTION 

The investigation of the topological and metric prop- 
erties of potential energy surface (PES), often referred to 
as "energy landscape" , is a useful and powerful tool for 
studying slow dynamics in condensed matter, especially 
in those cases where the lack of order (as for example 
in supercooled liquids) inhibits the use of the analytical 
tools pertaining to the crystalline state 0, 0, 0, 0, IE El • 
The PES approach has been successfully applied to the 
study of many different interacting systems (glasses, pro- 
teins, sheared materials, and so on). The PES approach 
started with the introduction of the fruitful concept of 
inherent structures 0. In the last years, several steps to- 
ward a more detailed description of the statistical prop- 
erties of the PES have been performed, most of them 
pointing toward a better understanding of the relation- 
ship between the landscape properties and the emergent 
dynamical behavior of the analyzed systems. 

Among others, two landscape-based approaches have 
proven to be particularly stimulating. The first one 
concerns the detailed analysis of the inherent structures 
(i.e. the configurations at the minima of potential en- 
ergy) visited by the system at different temperatures. 
This method has allowed to clarify many interesting phe- 
nomena, as, for example, the thermodynamic picture of 
the supercooled liquid regime based on the configura- 
tional entropy Q , the relationship between fragility and 
properties of inherent structures |5|, the analysis of dif- 
fusion processes in terms of visited inherent structures 
HHEl, or the interpretation of the effective fluctuation- 
dissipation temperature in the out-of-equilibrium re gim e 
in terms of inherent structures visited during aging , 
only to cite a few. The second approach is based on the 
analysis of the eigenvalues (normal modes) of the Hessian 
at the instantaneous configurations during the dynamic 
evolution of the system, from here the name instanta- 



neous normal mode approach (for an introduction and 
an extended application of this method see the works of 
Keyes and coworkers |0, El). This approach allowed 
to relate the emergent diffusive processes to the features 
of the landscape, opening the way to the interpretation 
of diffusion in terms of accessible paths in the multidi- 
mensional energy surface. Promising steps was obtained 
i) using simultaneously both the instantaneous normal 
mode approach and the inherent structure one, in order 
to identify the relevant slow diffusive directions 0, ^| , 
and ii) by analyzing the reaction paths in order to elim- 
inate the non-diffusive unstable modes |16| . 

Recently, a further approach has been introduced 17, 
[T3l and applied to the study of supercooled liquids 0,|2(J, 
HlH HI liH]. This approach is based on the anal- 
ysis of the saddles of the potential energy surface and 
has provided new insight in the analysis of the dynamic 
crossover taking place on lowering the temperature in 
supercooled liquids. Indeed it allows to characterize the 
dynamic transition temperature Tmct (mode-coupling 
temperature |26|) as the temperature where the order 
(fractional number of negative eigenvalues of the Hessian 
matrix) of the saddles vanishes. This finding suggested 
the following scenario for the dynamics: above Tmct 
the representative point in the configuration space lies 
close to the saddles and the relevant dynamic process is 
the diffusion among multidimensional saddle points, i.e. 
the diffusion takes place along paths at almost constant 
potential energy, and the limiting factors to particles dif- 
fusion are "entropic" , rather than "energetic" , barriers. 
On the contrary, below Tmct the minimum-to-minimum 
diffusion processes dominate and the "true" barrier jump 
controls the diffusive dynamics. A clear landscape-based 
interpretation of the dynamic behavior of the system is 
then provided. It is important to mention here that the 
term "saddles" is not mathematically correct, as the way 
the saddles have been defined in Ref.s [171 UM is based 
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on the partition of the configuration space in basins of 
attraction of the minima of the "pseudopotential" |2]j 
W= | VV^| 2 . It is clear that the absolute minima of W, lo- 
cated at W — 0, are true saddles of the energy surface (for 
simplicity of notation we call "saddles" also the minima 
and maxima of V) while the local minima of W (those 
with W > 0) correspond to points with (at least) one 
inflection direction, and are not saddles in mathemati- 
cal sense, rather they are "shoulders" along the inflec- 
tion direction. As pointed out by Doye and Wales [20| . 
the local minima of W , and not the absolute ones, are 
very often encountered during the minimization proce- 
dure. However, as it will be clear soon, the properties of 
the local and absolute minima of W which are actually 
important in determining the diffusive behavior are ex- 
actly the same. For this reason we call the local minima 
of W quasisaddles, to emphasize the fact that they carry 
the same information as saddles, even if they are geomet- 
rically different in nature (for a more detailed discussion 

see [2H E3)- 

Besides the landscape picture of the dynamic pro- 
cesses, the study of saddles has also permitted a quanti- 
tative characterization of the main features of the PES of 
liquid systems. Indeed, important PES properties, as the 
mean energy elevation of saddles from underling minima 
or the Euclidean distances among saddles, can be inferred 
from the analysis of saddle properties. It emerges an high 
regularity of the PES, with few parameters describing the 
spatial and energetic location of saddles. 

In this work we apply the saddle-approach to differ- 
ent model liquids (Lennard- Jones like pair potentials), in 
order to better understand the relationship between land- 
scape properties and slow dynamics, and in order to evi- 
dencing the existence of general features of the PES. The 
main result of this work is the existence of master curves 
both for temperature-dependent properties (saddle order 
vs. T) and for landscape properties (saddle energy vs. 
order), once energies and temperatures are normalized 
to Tmct- This is a very strict relationship between dy- 
namics and landscape features: differences in the PES 
for different systems simply define different Tmct val- 
ues, and once scaled by these values, one obtains exactly 
the same behavior. In other words, it appears that the 
PESs are very similar, the only differences being the val- 
ues of few parameters describing them (like the mean 
elevation barriers AE - mean elevation of saddles of or- 
der one from underlying minima) that lead to different 
values of dynamical quantities (Tmct)- The last point is 
of particular importance: for all the systems investigated 
we obtain that the value of Tmct is about 1/10 of the 
energy barrier AE, suggesting a kind of universality in 
the rearrangement processes governing the diffusion. 



II. MODELS 

We numerically investigated four different Lennard- 
Jones like model systems, all composed of iV = 256 parti- 



cles inside a cubic box with periodic boundary conditions. 
These are: 

1. the modified monatomic Lennard- Jones (ML J) 
|3lj. at ,0=1.0 (hereafter all the quantities will be 
expressed in LJ reduced units) 

V ML j(r) = 4e [(a/r) 12 - [a/rf] + SV , (1) 

where SV is a (small) many-body term that inhibits 
crystallization 

SV = aHf d(S(q) - So) [S(q) - So} 2 . (2) 

S(q) is the static structure factor, the sum is 
made over all q with q max —A<\q\<q max +A, where 
1max—7-12(p) 1 ^ 3 and A=0.34, and the values of the 
parameters are a=0.8 and So=lG. 

2. the modified monatomic soft spheres (MSS), at 
p=l.Q, 

V MSS (r) =4e(o-/r) 12 + SV , (3) 
where SV is defined in Eq. [3 

3. the binary mixture Lennard- Jones 80-20 (BMLJ) 
[12, at density ,0=1.2, 

Vbmlj{t) = 4e af3 [{(J a p/r) 12 - [a a p/rf] , (4) 

where the values of the parameters are those of the 
Kob- Andersen mixture (0x4 = 1, uab =0.8, obb = 
0.88, e AA = l, e A s = l-5, e BB = 0.5); 

4. a variant of the binary mixture Lennard- Jones 
(BMLJ 2 ), at p=1.2, in which the values of <jaa and 
obb were exchanged. 

In the case of BMLJ and BMLJ2, the interaction poten- 
tial is tapered at long distances between r± = IA'&oaa _ 
r < 2.56er aa = T2 with the following fifth-order smooth- 
ing function T[r) = 1 + (fi — r) 3 (6r 2 + (3r + ri)(ri — 
5r 2 ) + I0ri)/(r 2 - n) 5 . In this way the potential, the 
forces and their derivatives are continous, the energy can 
be kept constant to better than 1/10 5 over 100 millions 
of time steps. The MLJ and MSS potential have been 
simply cut and shifted at 2.5a. 

We performed standard molecular dynamics simula- 
tions at equilibrium (NVE ensemble), in a temperature 
range from T — 2 down to the lowest temperature that 
can be equilibrated in the MD run (this temperature 
is strongly model dependent). Along the equilibrium 
molecular dynamics trajectories at a given temperature 
we analyzed the properties of i) the instantaneous config- 
urations; ii) the inherent structures (minima); and iii) the 
saddle configurations. About 1000 configurations have 
been analyzed for each temperature and for each system. 
The inherent structures associated to instantaneous con- 
figurations are obtained by a conjugate-gradient mini- 
mization procedure on the total potential energy. For 
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Models 


P 


Tmct 


AE 


AE* 


MLJ 


1.0 


0.475 


4.43 


9.3 


MSS 


1.0 


0.210 


2.06 


9.8 


BMLJ 


1.2 


0.435 


4.16 


9.6 


BMLJ 2 


1.2 


0.605 


5.93 


9.8 



TABLE I: For the different Lennard- Jones like models we re- 
port the investigated density p, the mode-coupling tempera- 
ture Tmct (estimated from the apparent power-law vanish- 
ing of the diffusion coefficient), the mean barrier values AE 
(mean elevation of order-one saddles from underlying min- 
ima) and the reduced barrier height AE* — AE /Tmct ■ All 
the quantities are in LJ reduced units. 



saddles, a similar minimization procedure has been ap- 
plied to the pseudo-potential W = | W | 2 . The tapering of 
the BMLJ and BMLJ2 potentials allows the minimization 
procedures of both V and W to work correctly as they 
are not affected by small discontinuities in the derivative 
of V and W. The importance of avoiding discontinu- 
ities in order to obtain good W minimization has been 
recently underlined in Ref. |23j where the LBFGS algo- 
rithm [33| was used. However to obtain good minimiza- 
tions of W, even for a "small" system of 256 particles 
(i.e. 768 dimensions), is a stiff problem. We tested dif- 
ferent minimization algorithms (steepest-descent, Gauss- 
Newton, preconditioned conjugate gradient, Levenberg- 
Marquardt |3J|) but they eventually stick in some points 
of the configuration space, where the algorithm decrease 
more and more the step size, and the search becomes in- 
efficient and possibly stops. Different algorithms usually 
stick in different points. Sometimes the same algorithm 
who stuck in a given point can be effective in overcoming 
the critical situation if a larger step is used. Therefore, 
in the present work, a complex flow chart with various 
algorithms was used to obtain good minima (the details 
of the numerical algorithms will be presented elsewhere 
[35|). We want to remark that in this way the calls to the 
function W are always less than 3500 (average w 1500) 
and less than 1000 (average « 200) to the derivative of 
W. 

For all the analyzed configuration points (instanta- 
neous, minima and saddles) we store the energies per 
particle (e, e IS and e s respectively), and for instanta- 
neous and saddles we also determine their order n and 
n s , defined as the fractional number of negative eigenval- 
ues of the Hessian, i.e. the absolute number of negative 
curvatures over 3N (for inherent structures one obviously 
has n IS =0). 



III. SADDLES AND QUASISADDLES 

First of all we focus our attention on the differences 
between saddles (absolute minima of W) and quasisad- 
dles (local minima of W). As an example, Fig. n snows > 
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Log 10 (W/W mst ) 

FIG. 1: Histogram of the ratio W /Wi nst ,, i.e. the value at 
the minima of W with respect to the value at instantaneous 
configurations, at T = 2 for BMLJ2 (6000 configurations ana- 
lyzed). The higher region corresponds to quasisaddles (local 
minima of W), while the lower one to true saddles (absolute 
minima) . 



for the case of BMLJ2 model, the histogram of the value 
of the pseudopotential W at the minima (6000 config- 
urations analyzed at T = 2). The values of W at the 
minima are normalized to the values at the correspond- 
ing instantaneous configurations, i.e. to the value of W 
before starting the H^-minimization procedure. We ob- 
serve two very well distinct regions: the one with higher 
W values corresponds to local minima of W (quasisad- 
dles) , the lower one corresponds to absolute minima (true 
saddles). The non-zero values of W on the low-VF peak is 
due to the finite precision and/or threshold employed in 
the minimization procedure. A closer inspection of the 
eigenvalues of the Hessian shows that the quasisaddles 
are points with only one extra zero eigenvalue pfj| | (be- 
sides the three connected to the global translations), cor- 
responding to an inflection one-dimensional profile along 
the corresponding eigenvector. The fact that in the plot 
the two regions are well separated, allows to discriminate 
true and false saddles in a clear way. On the contrary, 
no clear separation has been found between saddles and 
quasisaddles from the analysis of the eigenvalues: due 
to the finite precision the found eigenvalues relative to 
the inflection points are different from zero of the same 
amount of the lowest frequency eigenvalues of real vibra- 
tional (or diffusive) modes. As it is evident from Fig. ^ 
true-saddles are very few and their number are found to 
decrease on lowering the temperature (e.g. for BMLJ2 in 
Fig. [2 about 5% at T = 2, and for BMLJ about 2% at 
T = 2 and less than 1% at T = 0.48). 

An interesting observation arises from the analysis of 
the behavior of the T-dependence of the number of nega- 
tive curvatures in the " true" saddles and in the quasisad- 
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dies separately. In Fig. [5] the saddle order is shown as 
a function of the temperature using only the true- (full 
symbols) and the quasi- (open symbols) saddles, for the 
cases of BMLJ (triangles) and BMLJ2 (squares) models 
(we note that in the BMLJ2 case, due to the appearance 
of crystallization, the data are available only for T > 1). 
The coincidence between the two set of data indicates 
that, as far as the temperature dependence of their char- 
acteristics (order and energy) is concerned, quasisaddles 
and true saddles share the same properties. Also other 
properties, as for example the spectral features (i.e. the 
density of vibrational states), of quasi-saddles and true- 
saddles are found indistinguishable |25j. This finding 
suggests that, no matter if saddles or quasisaddles, the 
minimization of W leads to points of the PES that are 
relevant for a landscape-based interpretation of the slow 
dynamics of the system: the order extrapolates to zero 
at the mode-coupling temperature Tmct (see Table [I] for 
the values of Tmct, estimated from diffusivity data, for 
the different models), indicating that at this tempera- 
ture the properties of the landscape probed by the sys- 
tem manifest a kind of discontinuity (the number of open 
directions, related to the saddle order, goes to zero and 
the dynamical processes change their characteristics) . In 
other words, the minimization of W seems to be a good 
method to get ride of the fast degrees of freedom and to 
keep information only on the slow degrees relevant for 
the slowing down of the dynamics taking place in super- 
cooled regime. 



IV. GENERAL FEATURES OF THE PES 

We now turn our attention to the existence of common 
features among the different model systems analyzed. 



A. T-dependent properties 

As already pointed out in Ref. [iH |28[. the (quasi-) 
saddle order, n s , vanishes as T approaches Tmct from 
above. At a first sight, it seems that the specific behav- 
ior of n s (T) is a model-dependent property (see Fig.|2J). 
However, we observe that after the scaling the tempera- 
ture scale by a specific sample dependent quantity, i. e. 
by Tmct, all the models behave similarly. In Fig. [3] the 
saddle order n s is reported as a function of reduced tem- 
perature T/Tmct- All the curves for the different sys- 
tems collapse into a single master curve. The latter can 
be fitted by a power law 

n s = n (—?- l) , (5) 

\J-MCT J 

with 7 = 0.85 and n — 0.025 (in the fitting procedure, the 
values of Tmct, reported in Tabled are kept fixed to the 
ones derived by the fit of the power-law behavior of the 
diffusion coefficient). A similar master plot is obtained 
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FIG. 2: Temperature dependence of the fractional order of 
true saddles (full symbols)and quasisaddles (open symbols), 
for BMLJ (triangles) and BMLJ2 (squares). Dashed lines are 
power law fits. 



also for the relation between the saddle energy and the 
temperature. These results suggest a universal behavior 
(at least for the LJ-like model systems analyzed here): 
at a given reduced temperature T* = T /Tmct all the 
systems visit saddles with the same properties (hereafter 
we will indicate with "*" the temperature and the energy 
scaled by Tmct)- One could conjecture that this univer- 
sality is due to the repulsive part of the pair potential 
r~ 12 (common to all the systems), that dominates over 
the attractive one at the studied densities. However, the 
facts that the curves superimpose each other quite well 
in the whole temperature range and that non-LJ systems 
(as, for example, the Morse potential - see the next sec- 
tion) show a similar behavior, seems to indicate that the 
observed universality is not trivially related to the repul- 
sive part of the interaction potential. Finally, we want 
to remark that the small value of n indicates that even 
at temperature twice that of the MCT critical point, the 
system is visiting saddles of low order (n s ~ 0.025) , so 
indicating that at T = 2T M ct the closest saddle, accord- 
ing to the partitioning defined by the minimization of W, 
is far below the top of the landscape. 



B. Energy barriers and Tmct 

The existence of common and general features of the 
PES emerges in a clear way from the comparative analy- 
sis of the energy and of the order of the saddles. In Fig.|U 
part A the energy elevation Ae s =e s — e IS of the saddles 
from the underling minima is plotted as a function of 
the saddle order n s for the different investigated models. 
As already observed 01 , there exists a proportionality 
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FIG. 3: Saddle order n s as a function of reduced temperature 
T/Tmct, for all the analyzed systems. The dashed line is a 
power law with exponent 7 = 0.85. For MLJ and MSS p— 1.0, 
while for BMLJ and BMLJ 2 p = 1.2. 



between these two quantities, indicating a simple organi- 
zation of the PES: saddles are equally spaced in energy 
over the minima. The slopes of the different straight lines 
in Fig. 0] determine the elementary energy elevation AE 
of saddles of order n from saddles of order n — 1: 



3 dn 



(6) 



where the factor 3 is due to the fact that energies are 
per particles (N) and the fractional order per degrees of 
freedom (37V). The values of AE obtained for the various 
systems are reported in Table [I] A possible explanation 
of the linear relationship observed in Eq. is that there 
exist in the system several spatially uncorrelated rear- 
ranging regions, each experiencing a mean barrier energy 
AE. In other words, if the system as a whole lies on a 
saddle of order to, this is due to the fact that there are 
to uncorrelated subsystems each one visiting a saddle of 
order 1. The analysis of the specific atomic motion as- 
sociated to these saddles, needed to assess or disprove 
the validity of this hypothesis, is beyond the aim of the 
present work. 

A very interesting and surprising result is obtained by 
scaling the energy values reported in Fig. 0] part A to 
the mode-coupling temperature Tmct (Kb — 1), obtain- 
ing again a single master curve (see Fig. 0] part B). The 
landscapes of different systems seem to share common 
features, with only one parameter describing the orga- 
nization of saddles, i. e. the mean elevation AE, that 
becomes an universal parameter (AE* = AE /Tmct — 
9 10) once normalized to the mode-coupling temper- 
ature (see the last column of Table QJ. In other words, 
all the models have the common property that the el- 
ementary barrier height is about 10 times the critical 
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FIG. 4: A) Energy elevation of saddles from underling min- 
ima e s — e IS against saddle order n 3 ; B) Energy elevation 
rescaled by mode-coupling temperature Tmct against saddle 
order n s . Dashed straight line is a guide to the eyes. 



temperature Tmct- 



AE ~ 10 Tmct . 



(7) 



This relation have been numerically proved for the four 
potential models investigated here. The same relation 
also holds for another LJ-like model, the binary mixture 
soft-sphere model (BMSS) investigated in Ref. [24| (at 
p = 1). This observation gives further support to the 
universality of Eq. If this is a particular characteristic 
of Lennard-Jones like models or a general feature of a 
more wide class of simple liquids is a open and interesting 
question which remains to be answered. 

We can try to give a first answer to this question an- 
alyzing the available data in the literature for other sys- 
tems. To our knowledge, besides the Lennard-Jones like 
systems, a saddle-based analysis has been performed only 
for the Morse potential ■ The Morse potential, used in 
Ref. H2, is denned as V a (r) = e[l-exp(a(l-r/r e ))] 2 -e, 
where e is the well depth, r e is the intcrparticle distance 
and the parameter a is inversely correlated to the range 
of the potential. Differently from soft-spheres and LJ, 
the Morse potential is finite as r — > 0. Unfortunately, 
equilibrium simulation based on the Morse potential are 
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difficult, since the undercooled system crystallizes eas- 
ily. Therefore simulations reported in Ref. [22j have been 
performed only well above Tmct- It was found that the 
larger the value of a is , the further the distance between 
the temperature of the lowest non-crystalline simulation 
and Tmct is. In this study, a linear dependence between 
n s and e s have been observed and the values of de s /dn s 
normalized to Tmct are in agreement with Eq. [7] for the 
three smaller a values a = 4, 5 and 6 (AE* are in the 
range 9.3-^10.5 H^|). For the two highest a values, a = 9 
and 12, the reported values for AE* are quite different 
(4.6 and 3.5 respectively). Further studies, for example 
focussing on binary mixture systems, are requested to 
find out if such discrepancy is due to an approximate 
determination of Tmct for a = 9 and 12 (which was ob- 
tained by extrapolating n s (T) from a temperature region 
where n s is far away from zero, n s (T) > 0.2). Uncer- 
tanties in the estimates of Tmct at large a are also con- 
sistent with the unexpected non-monotonic dependence 
of Tmct with a reported in Ref. [22^. We conclude that, 
for all a values for which the reliability of the data is un- 
questionable, the Morse potential landscape shares the 
same characteristic of those of the LJ-like potentials. 

In all the other model systems studied in the literature 
we do not have a direct information on the saddle energy 
elevation. However, the existence of a well defined bar- 
rier energy scale AE in the PES is expected to control 
the activation processes at low temperature, giving rise 
to an Arrhenius behavior of the transport properties at 
temperatures below Tmct- The existence of Arrhenius 
law in LJ-like systems - that are basically "fragile", in 
the Angell classification scheme Q - would be, per se, 
surprising (however, the degree of fragility of LJ systems 
is a matter of debate |37jl- 

The simulations below Tmct are very difficult to per- 
form, due to the extremely long relaxation times in this 
regime and a direct inspection of the expected Arrhenius 
behavior is not easy to pursuit. Only very recently such a 
kind of analysis has been performed for the BMLJ model 
at p = 1.2 38] . In that work an Arrhenius behavior was 
actually found in the temperature dependence of the dif- 
fusion coefficient below Tmct ■ T) oc exp(— AEatt/T) 
(we use the symbol AEatt for the activation energy in 
the Arrhenius law of the diffusivity, to distinguish it from 
the energy barrier AE determined from the saddles anal- 
ysis of the PES), with a value of AEatt — 8.1. The 
observed Arrhenius behavior is somewhat surprising in 
this "fragile" liquid models, and seems to indicate that 
close to Tmct activated processes start to be relevant 
and dominate the dynamics. However, the value of the 
activation energy AEatt found in Ref. |3S| is not equal 
to the elementary barrier energy AE estimated from the 
saddles analysis (see Table QJ, but it is about twice that 
value: AEArr/AE ~ 1.9. Re-analyzing our data for 
the MLJ model (for which we have few thermodynamic 
points equilibrated close to but below Tmct), we find 
that the above reported ratio is compatible with MLJ 
data (see Fig. |SJ), even if statistic is poorer than that 
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FIG. 5: Diffusivity D as a function of inverse temperature 
1/T for the MLJ model. Straight line represents the mode 
coupling like power law fit. Dashed line is the Arrhenius law 
with energy barrier AEatt — 1.9AE — 8.4 (AE is the energy 
barrier from saddles - see TableHJ , following the corresponding 
relation obtained for the BMLJ case (Ref. |3gln. 



of BMLJ case and the equilibrium condition is not fully 
satisfied by the lowest two temperature points. If such 
observation has general validity, then Arrhenius behav- 
ior should be observed below Tmct, with an activation 
energy value about 2 times the value of the elementary 
saddle energy barrier [so obtaining a value of reduced bar- 
rier energy (normalized to Tmct) A£^ rr ~ 18-^20]. The 
origin of this factor two needs to be further clarified. To 
this aim, it is important to underline that the " effective" 
energy barriers for activated processes as seen by the dy- 
namics (i. e. those entering in the Arrhenius law for the 
diffusivity) can be higher than the minimum-to-saddle 
energy difference (as measured directly by analyzing the 
PES). This can be due to the fact that the true diffusive 
path in the landscape 39] could pass higher in energy 
with respect to the saddle point, in order, for example, 
to minimize the minimum-to-minimum path length (i. e. 
for entropic reasons). In this respect, it is worth to men- 
tion that a non-coincidence between the relaxation times 
determined either through MD simulations or through 
the direct inspection of the PES has been observed in 
the simulation of a model protein during the folding pro- 
cess 0]. I n particular, the results in Ref. 0] indicates 
that the effective saddle height is larger than the actual 
one. 

Having in mind that AE* Arr ~ 2AE* and that AE* ~ 
10 (i.e. AE* Arr ~ 20), we can try to analyze what is 
observed for other model potentials existing in the liter- 
ature where the D(T) has been determined. We found 
three different models for which a low temperature anal- 
ysis of D(T) has been performed via molecular dynam- 
ics: i) The BKS-silica model 01 , f° r which the values 
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Models 


AE* 




MLJ 


9.3 


17.7 


MSS 


9.8 




BMLJ 


9.6 


18.6" 


BMLJ 2 


9.8 




BMSS 6 


9.1 




Morse c 


9.3 ~ 10.5 




Silica (BKS) d 




164-18 


OTP c 




20 4 28 


Water (SPC/E) / 




40 



°Ref. Fil 
6 Ref. pi 

c Ref. |2l| (obtained for a = 4, 5,6) 
d Ref. EJ 
e Ref. Eg 
^Ref. EJ 

TABLE II: Reduced energy barrier heights estimated from 
saddles (AE* = AE /Tmct) and from low-temperature Arrhe- 
nius law of diffusivity (AE Arr = AEatt/Tmct) for different 
model systems. The data of MLJ, MSS, BMLJ and BMLJ 2 
are from this work (except the AE\ rr for BMLJ, that is from 
Ref. H3). 

of AE* Arr are 16.2 and 18.0, for the self diffusion of O 
and Si respectively; ii) the Lewis and Wahnstrom ortho- 
terphenyl model |42|, for which the temperature depen- 
dence of the molecular center of mass diffusion coefficient 
at five different densities give values of reduced barrier 
energy AE* Arr ~ 20 4 28 (except the lowest density that 
gives a value of about 10); in) The SPC/E- water model 
|43|. for which one finds AE Arr ~ 40. Table QT| sum- 
marizes the known results on energy barrier heights esti- 
mated from saddles and from Arrhenius low-temperature 
dependence of diffusivity. The values for MLJ, MSS, 
BMLJ and BMLJ2 are from the present work, except 
the AE* Arr for BMLJ that is obtained from Ref. In 
future works we will try to determine the saddle-barriers 
AE* for non-LJ systems (the last three systems in the 
Table), in order to have a better understanding of the di- 
versity of the different landscapes. In conclusion, besides 
the case of water, the other systems seem to be in agree- 
ment with the findings of this work (the values of the 
reduced barrier energies are of the same order), evidenc- 
ing a quite general universality of the observed relations. 
A deeper understanding of the differences among various 
model liquids deserves further investigations. 



present work can be summarized as: 

• the coincidence between the temperature depen- 
dence of the quasisaddles and of the true saddles 
properties; 

• the existence of master curves for saddle proper- 
ties, once energies and temperatures are rescaled 
by the mode coupling critical temperature Tmct', 

• the existence of a universal relationship between 
the mode-coupling temperature and the mean 
energy barrier height AE ~ 10 Tmct, that seems 
to extend beyond the class of the Lennard- Jones 
like models analyzed here. 

Finally, we would like to point out that it already 
exists in the literature an hint on the existence of a linear 
relationship between AEa tt and the mode coupling 
critical temperature. Indeed, in a large class of glassy 
system one experimentally observes a linear relationship 
between the glass transition temperature T g and the 
infinite- frequency shear modulus Goo |H|: T g oc Goo- 
If we use the findings of our work (AE* Arr ~ 10, 
i.e. AEa„ ~ 10 Tmct) and we allow ourselves to 
confuse T g with Tmct, the following relation emerges: 
AEatt oc Goo- This relation is the prescription of the 
"shoving model" introduced thirty years ago by Nemilov 
|45| and recently put in a more rigorous form by Dyre 
et al. |46l |. The validity of the proportionality between 
AEArr and Goo has been proved for different glasses, 
and, together with the linear relationship between T g 
and Goo, give further support to the finding of the 
present work, i. e. the apparent universality of the ratio 

AEatt/TmCT- 
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V. CONCLUSIONS 



In conclusion, despite complex and disordered in na- 
ture, the simple liquid PES seems to exhibit few general 
and regular features, useful both to bring important in- 
sight for the understanding of the relevant diffusion pro- 
cesses taking place in supercooled liquids and to con- 
struct simplified PES models. The main findings of the 
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